On the nonlinear stage of Modulation Instability 
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We study the nonlinear stage of the modulation instability of a condensate in the framework 
of the focusing Nonlinear Schrodinger Equation. We find a general A''-sohtonic solution of the 
focusing NLSE in the presence of a condensate by using the dressing method. We separate a special 
designated class of "regular solitonic solutions" that do not disturb phases of the condensate at 
infinity by coordinate. All regular solitonic solutions can be treated as localized perturbations of 
the condensate. We find an important class of "superregular solitonic solutions" which are small 
perturbations at certain a moment of time. They describe the nonlinear stage of the modulation 
instability of the condensate. 
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-Introduction - The focusing Nonlinear Schrodinger 
Equation (NLSE) is a universal model for studying quasi- 
monochromatic wave propagation in weakly nonlinear 
media. In particular it describes waves on deep water 
1 1 , waves in optical fibers Q , Langmuir waves in plasma 
3| . Also this equation is a model of weakly interacting 
Bose gas with attraction between particles. 

The NLSE has a simple solution, the monochromatic 
wave with frequency depended on amplitude - the con- 
densate. The condensate is unstable with respect to mod- 
ulation instability (The history of modulation instability 
is described in What is a nonlinear stage of mod- 
ulation instability? In spatial dimension D — 2,3, the 
answer is known - modulation instability leads to forma- 
tion of finite time singularities - collapses Q . 

In dimension D = 1 collapses are forbidden. However 
in this case development of modulation instability leads 
to formation of "extreme" (rogue, freak) waves where 
energy density exceeds the mean level by order of magni- 
tude. As a result, the study of long-time consequences of 
modulation instability is a problem of big practical im- 
portance, crucial for creation of a rogue wave theory in 
the ocean and a theory of extreme events in optical lines. 

NLSE is a completely integrable system having 
many exact solutions. It is natural to hope that the 
nonlinear development of the modulation instability is 
described by such a solution. 

In what follows we speak only about instability grow- 
ing from a small spatially localized perturbation of con- 
densate. Historically the first such solution was found by 
Peregrine in 1983 [yj. This equation attracted a lot of at- 
tention Its experimental confirmation was claimed 
d, 0]. In 1985 a second order Peregrine solution was 
found [l3] ■ Today " multi-Peregrine" are actively studied 

All these 
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by different groups (see for instance 
solutions have a weak point - they are small perturba- 
tions of the condensate only in the limit i — > — oo. These 
solutions are homoclinic - they describe freak waves ap- 



pearing "from nowhere" and completely disappearing in 
the future. Meanwhile numerical modeling of modula- 
tion instability 13] as well as numerical solutions of ex- 
act Euler equation for potential flow of ideal deep fluid 



14j demonstrate formation of propagating oscillating so- 
lutions ("breathers"). Thus the "homoclinic scenario" of 
modulation instability development can be disputed. 

In this paper we present another class of exact solutions 
of NLSE which are small perturbations of condensate not 
at t ±oo but in the initial moment of time t = 0. This 
is the special class of 2A^-solitonic solutions of NLSE in 
the presence of a condensate. 

A solitonic solution of the NLSE in presence of conden- 
sate was first found by Kuznetsov in 1977 15] and redis- 
covered later in Il7| . An important solitonic solution 
was found in 1985 by Akhmediev, Eleonskii and Kulagin 
[13 (the so called " Akhmediev breather" ) . This is a solu- 
tion periodic in space and homoclinic in time. The Pere- 
grine solution is the limiting case of both Kuznetsov's 
and Akhmediev's solutions. More general solitonic solu- 



tions were found in |18l - l22l | . To construct exact solutions 
of NLSE describing the development of a localized small 
perturbation of the condensate one has use 2iV-solitonic 
solutions. There are several mathematical schemes for 
constructing such solutions (Hirota method, method of 
Darboux transformation etc.) We prefer the "dressing 



method" elaborated in [23|. The details of our mathe- 
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matical procedure can be found in 

The main result of our parer is the following. One can 
construct a broad class of 2N - solitonic solutions de- 
pending on 7N parameters which form at t = a small 
perturbation of condensate. Evolution of this perturba- 
tion lead to the formation of complicated "integrable tur- 



bulence" [25| where local concentration of energy easily 



exceeds in order of magnitude the energy density in the 
condensate. In the presented paper we describe only the 
simplest 2-solitonic solution with 7 free parameters which 
is at i = a small perturbation of condensate. 
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-NLSE- We write the NLSE in the following form 
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(1) 



The trivial condensate solution of equation ^ \s ip = 
ipo — A. One can consider A to be real. This solu- 
tion is unstable with respect to small perturbations. The 
growth rate of instability is T{p) = ^\f^A^^~^y^ ^ where 
•p is wave number of perturbation. In what follows we 
use the NLSE with non-vanishing boundary conditions 
— !> \A\^ at X — >■ ±cxD. Equation ([T]) is the compati- 
bility condition for the following_overdetermined linear 
system for a matrix function 4* 



Here 
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Here A is a spectral parameter. Suppose we know a 
certain particular solution ipo of equation ([I} together 
with the fundamental matrix solution 4'o(x,t, A) of sys- 
tem ([S]),©. Then one can construct a new solution 
If of equation ([T]) using the following recipe. Choose 
N complex numbers Afc (k=l,..N) ReXk > and an- 
other set of arbitrary complex numbers Ci, ..Cn- Denote 
Fk = ^o(a;,i, —^k) define N vectors g„ by relation 



Then a new solution is given by expression 
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Here M^p (a = 1, 2) is the following determinant 
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a/3 



qi.p 



Mi 
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Where Mnm is a Hermitian matrix: 



An + A* 



Mention that transformation g„ dnlm where a„ are 
arbitrary complex numbers, does not change the result 
of dressing. 

lffo = A 



Here 



s(A) -e- 



{x,t,\) 



(10) 



^kx + nt, k^^X^-A^, n = -iXk, s = - 



A 



X + k 



^0 has a cut from —A to A. We perform the Jukowsky 
transform and map this plane onto the outer part of the 
circle of unit radius. 



A = |(e + ^) e + r^o 



and use parametrization 



After redefinition of phase factor 
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Here 



w„ 



IV„ 



Vn = knX - iOnt - -On (14) 



and 



ffi„ — Asinh(z„) cos(q!„) 
kn — ^ cosh(z„) sin(Q;„) 

U = --^ cosh(2z„) sin(2a„) 
^2 

w„ = — sinh(2z„) cos(2a„) 



(15) 



If n = 1 we get a one-solitonic solution characterized 
by four parameters R = > l,a,9,fj,. The first two 
parameters define the location of the complex spectral 
parameter A which is actually a pole of the "dressing 
function" (see [2J]). In the Kuznetsov case a = 0, R > 1 
and the pole is located on the real axis outside of the cut. 

In the Akhmediev case the pole is located on the unit 
circle R — l,a 0. For the Peregrine solution R = 
l,a — 0. Now the pole is posed exactly in the branch 
point. In a general case i? > 1, a ^ 0, the pole is located 
on the complex plan outside the unit circle. Parameters 
9, jjL define the location and phase of the soliton. If 9 ~ 
0,/x — the soliton satisfies the symmetry condition. 



M = det{Mnm) 



(9) 



(p{-x,-t) = ip*ix,t) 



(16) 
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In the general case the sohton is filled with carrying wave 
moving with phase velocity. 



Vph = — - 



(17) 



The soliton's envelope moves with the group velocity 
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(18) 



In the Kuznetsov-Ma and Peregrine case Vph = 0. In the 
Akhmediev case Vph = 0, Vgr — oo. We are specifically 
interested in " near- Akhmediev soliton" when the pole is 
close to the unit circle (z << 1). This is a quasi-periodic 
solution of large size. 

L ~ ^ , \ , (19) 

2zAcos{a) 

It moves with high group and small phase velocity 
A 2Azcos(2q!) 



V„ 



Vph 



y group ~ • / \ • pn — • / \ 

sm(a)z sm(a) 
The "near- Akhmediev" solution is plotted on Fig. [TJ 

1^1 J,A=1 
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FIG. 1: The " near- Akhmediev" solution. (7?= 1.02, q = |) 

Note that the one-solitonic solution is defined by one 
vector 1^ ~ {qi,q2). Its shape is given by the following 
formula 



4cosh(2;)cos(a)("^| • "^2) 

^^^-^ itrP + lt.P 

This solution has the following asymptotics 

A general A''-solitonic solution has asymptotics 
(/J Ae^'^^'^'^i X ^ ±00 
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If we are interested in A^-solitonic solutions localized in 
a finite domain of space and not perturbing the remote 
condensate we must put 
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We call this solution a "regular solitonic solution". If we 
assume that the modulation instability develops from lo- 
calized perturbation, only a regular solution can be used 
as a model for its nonlinear behavior. 

Among one-solitonic solutions only the Kuznetsov-Ma 
and the Peregrine solutions are regular. In two-solitonic 
case we can construct a broad class of regular solutions. 

-Two-solitonic solution- The two solitonic solutions is 
defined by two vectors "^1 = (gii,gi2),Y2 = (921,922)- 
A general regular solitonic solution is depends on seven 
parameters Ri = e^^ i?2 = e^^, a, /.ti, /i2, 6*1, ^2- When 
zi — Z2 = z poles located in complex conjugated points. 
In this case twosolitonic solution can be presented in the 
following form. 



(f = A- 
sinh(z) sin(Q:)i?i 



^sinh(2z) sin(2Q;) x 
- i cosh(z) cos(Q;)i?2 



cosh (z) cos2(a)|S'|2 -|- sinh (z) sin {a)P 



(25) 



5 = (711522 - '7i2'72i P=\qi\\~¥2\ 

Bi = l"^il^(?2i'722 + |'^2pgngi2 

B2 ■ "^2)921 912 - {~ti ■ ~t 2)111122 



(26) 



Suppose that 5 7^ and z — 0. Then ip = A, thus the 
condensate is not perturbed. We call this phenomenon 
annihilation of Akhmediev breathers. If z — the an- 
nihilation is " not complete" . Solution is a small per- 
turbation at the moment of solitons collision. Let us 
put fJ-i ^ fJ-2 — 0. In this case solitons collide at 
[x — 0,t ~ 0). Let i?~l+e,z~£, then 



if = A + Sip 



5ip 



, M - iK 
-2eA—^ (27) 



M — 2 sin(a) sinh(C) sin(i9) cos ( — \ 
K = 2cos(a) cosh(C) cos(i9) sin^^^ + sin(0+) 
G = cos2 (^-^\ sin\?{C) + sin^ (^-^\ cosh^{C) 



d = 2Asin(Q!)a; — 

9+^01 + 



( = 2Ae cos(Q!)a; 

e- ^ei- 02 (28) 



li 6^ 7^ 0, that G and 6ip is a small perturbation of 
condensate. An initially small localized perturbation of 
the condensate generates a pair of near- Akhmediev soli- 
tons propagating with very fast group velocity in opposite 
directions. Typical development of these small localized 
perturbation of the condensate is presented on Fig. [2] 

Annihilation of solitons takes place for much more gen- 
eral class of solitonic solutions. Let us consider a 2A^ soli- 
tonic solution consisting of N pairs of Akhmediev solitons 
such that 



Rn = R 



■n+N 



1 



a 



-OLn+N 



(29) 
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t=0 



\<P\ ',A=1 
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cation to creation of statistical theory of rogue waves in 
ocean and optic fibers. 
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FIG. 2: Typical development of small localized perturbation 
of the condensate, a — l,e — 0.2, a — ^,61 = ^,62 



Here n = 1,...N. Each pair is characterized by addi- 
tional parameters jin, IJ^n+N ^(^mQn+N ■ We assert that if 
all 9^ = 6n + djy+n 7^ this solution annihilates com- 
pletely. If the annihilation is "incomplete" such a soli- 
tonic solution presents a small perturbation of conden- 
sate. We call this solutions a "superregular solitonic so- 
lutions" . Notice that an incomplete annihilation is not 
necessary symmetric. Even in the case = 1 one can 
put i?i = 1 -I- oie, i?2 = 1 + a2S, ai 7^ a2, e — !> 0. Such 
an initially small perturbation generates a pair of dif- 
ferent near Akhmediev solitons propagating in opposite 
directions. A generic solution of the mentioned type is 
a nonlinear superposition of N nonsymmetric pairs of 
near- Akhmediev solitons, and can be treated as a sort 
of " integrable turbulence" appearing as a result of non- 
linear development of the modulation instability. What 
is interesting is that when e — ;> this superposition is 
linear, so that at t = small perturbations generated 
by separate "superregular" pairs of soliton form an A^ - 
dimensional linear space. This remarkable fact will be 
discussed in another article. 

-Conclusion - We show that a broad class of small 
perturbations of the condensate leading to the develop- 
ment of modulation instability can be described by exact 
solitonic solutions. Our next plan is the statistical de- 
scription of complex solitonic solutions that appear as a 
result of modulation instability development with appli- 
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